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We provide the most precise measurement of the WW production cross section in pp coUisions 

to date at a center of mass energy of 1.96 TeV, and set limits on the associated trilincar gauge 
couplings. The WW Ivt'v {£,£' = e,/i) decay channels are analyzed in 1 fb""*^ of data collected 
by the DO detector at the Fermilab Tevatron Collider. The measured cross section is ij{pp 
WW) = 11.5 ± 2.1 (stat + syst) ± 0.7 (lumi) pb. One- and two-dimensional 95% C.L. limits on 
trilinear gauge couplings are provided. 



The non-Abclian gauge group structure of the eleetro- 
weak sector of the standard model (SM) predicts specific 
interactions between the 7, W, and Z bosons. Two ver- 
tices, WW-y and WWZ, provide important contributions 



to the pp — ^ WW production cross section. Understand- 
ing this process is imperative because it is an irreducible 
background to the most sensitive discovery channel for 
the Higgs boson at the Tevatron, H — > WW. A de- 
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tailed study of WW production also probes the triple 
gauge-boson couplings (TGCs), which are sensitive to 
low-energy manifestations of new physics from a higher 
mass scale, and is sensitive to the production and decay 
of new particles, such as the Higgs boson ilj. Study- 
ing WW production at the Fcrmilab Tevatron Collider 
provides an opportunity to explore constituent center 
of mass energies (VS) higher than that available at the 
CERN e+e- Colhder (LEP) 0], since SM WW produc- 
tion at the Tevatron has an average V§ — 245 GeV and 
a 57% probability for VS > 208 GeV The Tevatron 
experiments have been active in studying the WW cross 
section and TGCs in the past [1, 0, Q . In this Letter we 
present the most precise measurement of the WW pro- 
duction cross section in pp collisions to date and updated 
limits on anomalous WW^ and WW Z couplings. 

We examine WW production via the process pp — > 
W+W- l+vl'-D £' = e, /i; allowing for W tv ^ 
£ + 711/ decays) and use charged lepton transverse momen- 
tum (pt) distributions to study the TGCs. The decay of 
two W bosons into electrons or muons results in a pair 
of isolated, high-pr, oppositely charged leptons and a 
large amount of missing transverse energy (J^t) due to 
the escaping neutrinos. This analysis uses pp collisions 
at a center of mass energy of 1.96 TeV, as recorded by 
the DO detector Q at the Tevatron. A combination of 
single-electron (ee and e/j, channels) or single-muon (/i/i 
channel) triggers were used to collect the data, which 
correspond to integrated luminosities of 1104 (ee), 1072 
(e^), and 1002 (nfi) pb"! 0. 

Electrons are identified in the calorimeter by their elec- 
tromagnetic showers, which must occur within \ri\ < 1.1 
or 1.5 < \ri\ < 3.0 [8]. In the ee channel, at least one 
electron must satisfy |ry| < 1.1. Electron candidates must 
be spatially matched to a track from the central track- 
ing system, isolated from other energetic particles, and 
have a shape consistent with that of an electromagnetic 
shower. Electron candidates must also satisfy a tight re- 
quirement on a multivariate electron discriminant which 
takes into account track quality, shower shape, calorime- 
ter and track isolation, and E/p, where E is the calorime- 
ter cluster energy and p is the track momentum. The pr 
measurement of an electron is based on calorimeter en- 
ergy information and track position. 

Muons are reconstructed within |?7| < 2.0, must be 
spatially matched to a track from the central tracking 
system, and are required to have matched sets of wire 
and scintillator hits before and after the muon toroid. 
The detector support structure limits the muon system 
coverage in the region |?7| < 1.1 and 4.25 < 4> < 5.15 [1]; 
in this region a single set of matched wire and scintillator 
hits is required. Additionally, muons must be isolated 
such that the pt sum of other tracks in a cone TZ — 
yJ{l\'qY + (A0)2 < 0.5 is < 2.5 GcV and calorimeter 
energy within 0.1 < 7?. < 0.4 is < 2.5 GeV. 

The is determined based on the calorimeter energy 
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FIG. 1: Distributions of the (a) leading and (b) trailing 
lepton Pt after final selection, combined for all channels 
(ee -I- e/i 4- /i/i). Data are compared to estimated signal, 
CT(iyW) = 12 pb, and background sum. 



deposition distribution with respect to the interaction 
vertex. It is corrected for the electromagnetic or jet en- 
ergy scale, as appropriate, and the pt of muons. 

Signal acceptances and background processes are stud- 
ied with a detailed Monte Carlo (MC) simulation based 
on PYTHIA [§] in conjunction with the CTEQ6l1 [i3| par- 
ton distribution functions, with detector simulation car- 
ried out by GEANT [ll[. The Z boson pT spectrum in 
Z/7* ii MC events is adjusted to match data [l3 |. 

For each final state, we require the highest px (leading) 
lepton to have pr > 25 GeV, the trailing lepton to have 
Pt > 15 GeV, and the leptons to be of opposite charge. 
Both charged leptons are required to originate from the 
same vertex. The leptons must also have a minimum 
separation in rj-cj) space of TZee > 0.8 in the ee channel 
or Tlep,/p,^ > 0.5 in the efi and /i/i channels, in order to 
prevent overlap of the lepton isolation cones. 

Background contributions to WW production from 
M^-l-jets and multijet production are estimated from the 
data. Those from Z/7* U, tt, WZ, W-f, and ZZ are 
estimated from the MC simulation. 

After the initial event selection, the dominant back- 
ground in each channel is Z/j* ^ ££ {£ = e, n, t). Much 
of this background is removed by requiring > 45 (ee) , 
20 (e/i), or 35 (/i/i) GeV. For the ee channel, we require 
fr > 50 GeV if |Mz - mee| < 6 GeV to further re- 
duce the Z/'j* ££ background. In events contain- 
ing muons, a requirement on the azimuthal separation 
(A0) between the leptons is more effective at reducing 
the Z/j* — > ££ background than an invariant mass re- 
quirement, since the momentum resolution for high pT 
muons is poorer than the calorimeter energy resolution 
for electrons. The e/i channel additionally requires 
> 40 (instead of 20) GeV if A^g^ > 2.8, and the /i/i 
channel requires A0^^ < 2.45. 

Mismeasurement of the muon momentum can lead to 
spurious which is collinear with the muon direction. 
Especially in the /x/i channel, mismeasurement of the 
muon momentum can allow Z boson events to satisfy 
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the fx requirement. To suppress these events in the /i/i 
channel, we require that the track for each muon candi- 
date include at least one silicon microstrip tracker hit, 
for better momentum resolution, and that the azimuthal 
angle between each muon and the direction of the fx 
satisfies | cos(A0^^,^)| < 0.98. 

A second background is tt production followed by the 
leptonic decay of W bosons. This background can be 

suppressed by requiring qt ~ {pre + pre' +-^t| < 20 
(ee), 25 (e^), or 16 {^fi) GeV. This quantity is the pt of 
the WW system and is expected to be small for signal 
events. However, for ti production and other background 
processes, qr can be large, so this variable is a powerful 
discriminant against these backgrounds. 

The process is a background for only the ee and 
e^ channels, since the probability for a photon to be 
misidentified as a muon is negligible. We determine the 
probability that a photon is misidentified as an electron 
with photons from Z/j* eej decays and use it to 
correct the MC-based prediction of the W'y background. 
The H^+jets background, in which a jet is misidentified 
as an electron or muon, is determined from the data by 
selecting dilepton samples with loose and tight lepton re- 
quirements and setting up a system of linear equations 
to solve for the W+jet backgrounds after all event se- 
lection cuts, similar to the multijet background estima- 
tion performed in [l3l |. The multijet background con- 
tains jets that are misidentified as the two lepton candi- 
dates. It is represented by a data sample where the re- 
constructed leptons fail the lepton quality requirements. 
This sample is normalized with a factor determined at 
preselection using like-charged lepton events. It is as- 
sumed that misidentified jets result in randomly assigned 
charge signs. 

The leptonic decay of WZ and ZZ events can mimic 
the WW signal when one or more of the charged leptons 
is not reconstructed and instead contributes to.^^. The 
ZZ — > iiviy process is suppressed by the \Mz — rrieel or 
A(j)w cut. 

For each channel, the exact selection requirements on 
^T, Qt, and \Mz — irieel O'" ^i'U' chosen by per- 
forming a grid search on signal MC and expected back- 
ground, minimizing the combined statistical and system- 
atic uncertainty on the expected cross section measure- 
ment. The final lepton pT distributions are shown in 

Fig.[T][M]. 

The overall detection efficiency for signal events is de- 
termined using MC with full detector, trigger, and re- 
construction simulation and is 7.18% (ee), 13.43% (e^), 
and 5.34% (^^) for WW tvl'v {1,1' = e,/i) de- 
cays and 2.24% (ee), 4.36% (e^), and 1.30% {^i^i) for 
WW — > Tvlv Itvtv W + nv decays. The numbers of 
estimated signal and background events and the number 
of observed events for each channel after the final event 
selection are summarized in Table ID The observed events 



TABLE I: Numbers of signal and background events expected 
and number of events observed after the final event selection 
in each channel. Negligible contributions are not shown. Un- 
certainties include contributions from statistics and lepton se- 
lection efficiencies. 



Process 


ee 


ejj, 


MM 




0.27 ±0.20 


2.52 ±0.56 


0.76 ±0.36 


Z/^* TT 


0.26 ±0.05 


3.67 ±0.46 


— 


ti 


1.10±0.10 


3.79±0.17 


0.22 ±0.04 


W Z 


1.42 ± 0.14 


1.29 ± 0.14 


0.97±0.11 


ZZ 


1.70 ±0.04 


0.09 ±0.01 


0.84 ±0.03 


Wj 


0.23 ±0.16 


5.21 ±2.97 




W + iet 


6.09 ±1.72 


7.50 ±1.83 


0.12 ±0.24 


Multijet 


0.01 ±0.01 


0.14±0.13 




WW -* £1' 


10.98 ±0.59 


39.25 ±0.81 


7.18 ±0.34 


WW -+ eT/TT - 


1.40 ±0.20 


5.18 ±0.29 


0.71 ±0.10 


Total expected 


23.46 ± 1.90 


68.64 ±3.88 


10.79 ±0.58 


Data 


22 


64 


14 



are statistically consistent with the SM expectation in 
each channel. Assuming the W boson and t branching 
ratios from [l5| . the observations in data correspond to 
a(j)p^ WW) = 10.6±4.6 (stat)±1.9 (syst)±0.7 (lumi) 
pb in the ee channel, 10.8 ± 2.2 ± 1.1 ± 0.7 pb in the e/i 
channel, and 16.9 ± 5.7 ± 1.4 ± 1.0 pb in the fifj, chan- 
nel. The dominant sources of systematic uncertainty for 
each channel are the statistics associated with the esti- 
mation of the ly+jet contribution in the ee channel, the 
photon misidentification probability used to estimate the 
Wj contribution in the e/Lt channel, and the MC statistics 
for backgrounds in the /x/z channel [l3 | . 

The cross section measurements in the individual chan- 
nels are combined usin g th e best linear unbiased esti- 
mator (BLUE) method [Idl yielding: a{pp WW) = 
11.5 ± 2.1 (stat + syst) ± 0.7 (lumi) pb. The standard 
model calculation of the WW production cross section at 
the Tevatron center of mass energy is 12.0 ± 0.7 pb [13] ■ 

The TGCs that govern WW production can be pa- 
rameterized by a general Lorentz-invariant Lagrangian 
with fourteen independent complex coupling parameters, 
seven each for the WW^ and WWZ vertices [H . Limits 
on the anomalous couplings are often obtained by tak- 
ing the parameters to be real, enforcing electromagnetic 
gauge invariance, and assuming charge conjugation and 
parity invariance, reducing the number of independent 
couplings to five: ^f, kz, k^, Xz, and A-y (using no- 
tation from [1]). In the SM, gf — kz — — 1 and 
Xz = X^ = 0- The couplings that are non-zero in the SM 
are often expressed in terms of their deviation from the 
SM values, e.g. Agf = gf-l. Enforcing S'C/(2)i(^J7(l)y 
symmetry introduces two relationships between the re- 
maining parameters: kz = gf — {k-^ — l)tan^6vv and 
Xz = Xj, reducing the number of free parameters to 
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three [Ig. Alternatively, enforcing equality between the 
WW-i and WWZ vertices {WW-i=WWZ) such that 
— Kz, = Xz, and gi — 1 reduces the number 
of free parameters to two. 

One effect of introducing anomalous coupling parame- 
ters into the SM Lagrangian is an enhancement of the 
cross section for the qq Zj^* W^W~ process, 
which leads to unphysically large cross sections at high 
energy. Therefore, the anomalous couplings must vanish 
as the partonic center of mass energy \fh — > cxd. This is 
achieved by introducing a dipole form factor for an arbi- 
trary couphng a {gf , Ky, or Ay): a{s) ~ q;o/(1 + s/A^)^, 
where the form factor scale A is set by new physics, and 
limits are set in terms of Uq. Unitarity constraints pro- 
vide an upper limit for each coupling that is dependent 
on the choice of A. For this analysis we use A = 2 TeV, 
the approximate center of mass energy of the Tevatron. 

The leading order MC event generator by Hagiwara, 
Woodside, and Zeppenfeld [H is used to predict the 
changes in WW production cross section and kinematics 
as coupling parameters are varied about their SM values. 
At each point on a grid in TGC parameter space, events 
are generated and passed through a parameterized sim- 
ulation of the DO detector that is tuned to data. To en- 
hance the sensitivity to anomalous couplings, events are 
sorted by lepton pT into a two-dimensional histogram, 
using leading and trailing lepton pT values in the ee and 
/i/i channels, and e and /i pt values in the epL channel. For 
each bin in lepton pT space, the expected number of WW 
events produced is parameterized by a quadratic func- 
tion in three-dimensional (AK^,A-y,A(7f ) space or two- 
dimensional (Ak,A) space, as appropriate for the TGC re- 
lationship scenario under study. In the three-dimensional 
case, coupling parameters are investigated in pairs, with 
the third parameter fixed to the SM value. A likelihood 
surface is generated by considering all channels simul- 
taneously, integrating over the signal, background, and 
luminosity uncertainties with Gaussian distributions us- 
ing the same methodology as that used in previous stud- 
ies 0. 

The one-dimensional 95% C.L. limits for A = 2 TcV 
are determined to be -0.54 < Ak-^ < 0.83, -0.14 < 
A^ = Az < 0.18, and -0.14 < A^f < 0.30 under the 
SU{2)l ® [/(l)y-conserving constraints, and —0.12 < 
Ak^ = Ak^ < 0.35, with the same A limits as above, 
under the WWj=WWZ constraints. One- and two- 
dimensional 95% C.L. limits are shown in Fig.O 

In summary, we have made the most precise measure- 
ment of WW production at a hadronic collider to date, 
a{pp WW) = 11.5 ±2.1 (stat + syst) ±0.7 (lumi) pb, 
using 1 fb^^ of data at the DO experiment. This result 
is consistent with the SM prediction and previous Teva- 
tron results The selected event kinematics are 
used to significantly improve previous limits on anoma- 
lous TGCs from WW production at the Tevatron, re- 
ducing the allowed 95% C.L. interval for A^ = Xz and 
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FIG. 2: One and two-dimensional 95% C.L. limits when en- 
forcing SU{2) L ® U{X)y symmetry at A = 2 TeV, for (a) Ak^ 
vs. A-Y, (b) Ak^ vs. Agf , and (c) vs. Agf , each when 
the third free coupling is set to its SM value; limits when 
enforcing the WW'y=WW Z constraints are shown in (d). 
The curve represents the two-dimensional 95% C.L. contour 
and the ticks along the axes represent the one-dimensional 
95% C.L. limits. An asterisk (^) marks the point with the 
highest likelihood in the two-dimensional plane. 



Ak-^ — Akz by nearly a factor of two (Bl. [20|. 
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SUPPLEMENTAL MATERIAL 

The final Icpton pT distributions for individual analy- 
sis channels are shown in Figure [3l The uncertainty from 
each systematic source as a percentage of the final back- 
ground estimate for each channel is provided in Table [ill 

In Table |TT1 the "MC Statistics" uncertainty is based 
on the number of MC events used to estimate the num- 
ber of background events after final selection. The "MC 
Cross Section" uncertainty is based on the PDF and scale 
uncertainties associated with the theoretical cross sec- 
tions used to scale the MC-driven background estima- 
tions. The "MC Corrections" uncertainty accounts for 
data-driven reweighting of the MC that corrects for gen- 
eral event characteristics such as the primary vertex z- 
position and instantaneous luminosity distributions. The 
"7 e Mis-ID Rate" uncertainty is driven by statistics 
in the data sample used to determine that misidentifi- 
cation rate. Sources of systematic uncertainty for the 
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TABLE II: Uncertainty from each systematic source as a per- 
centage of the final background estimate for each channel. 



Systematic Source 






ee 


e/i 




MC Statistics 






2.6% 


3.6% 


12.8% 


MC Cross Section 






2.6% 


2.2% 


2.5% 


MC Corrections 






2.0% 


6.1% 


2.0% 


Electron ID 






0.3% 


0.3% 




Muon ID 








0.4% 


3.0% 


7 ^ e Mis-ID Rate 






1.2% 


12.1% 




W-l-jet Statistics 






15.5% 


7.5% 


6.8% 


W-l-jet e Efficiency <§. 


L Mis-ID 


Rate 


3.4% 


2.9% 




W^-l-jet /I Efficiency ^ 


z Mis-ID 


Rate 




1.4% 


2.1% 


Total Background Uncertainty 


16.4% 


16.4% 


15.3% 



M^-|-jet background estimation are separated into cate- 
gories based upon their correlations across analysis chan- 
nels. The "M^-l-jet Statistics" systematic is based on the 
number of events in data that pass a series of lepton iden- 
tification requirements in each analysis channel, indepen- 
dently. The 'W+iel Lepton Efficiency & Mis-ID Rate" 
uncertainties are based on measuring the rate at which a 
lepton passes loose and tight identification in Z/7* U 
candidate events in data and the rate at which a jet is 
misidentified as a lepton with loose or tight identification 
based on dijet data. 

The statistical uncertainty due to the number of ob- 
served events in each channel is the dominant source of 
final uncertainty. 
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FIG. 3: Distributions of the (a) leading and (b) trailing elec- 
tron pt in the ee channel, (c) eiectron and (d) muon pr in 
the e^l channel, and (e) leading and (f) trailing muon pr in 
the /i/i channel after final selection. Data are compared to 
estimated signal, (t(VFVF) = 12 pb, and background sum. 



